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■^f ■ 1. Introduction 

The problem of a one-dimensional harmonic oscillator is one of the important paradigms 
Q\ \ of the theoretical physics. Its solution in the classical approach is unique and simple, 
leading to the enormous applications in a wide range of the modern physics [1]. This 
problem in the non-relativistic quantum approach has at least the same importance, 
and probably the main reason is that it has very elegant solutions to the Schrodinger 
equation in terms of the well-known Hermite polynomials and the algebra of this problem 
can be easily factorized as Heisenberg algebra with its generators being the starting 
point of the quantum field theory. Here, one needs to note that there is another 
important solution to the Schrodinger equation in terms of the Hermite polynomials and 
it is the non-relativistic harmonic oscillator in an external homogeneous field [2]. The 
importance of this non-relativistic problem can also be observed through its numerous 
applications. For example, the theory of the Brownian motion of a quantum oscillator 
is developed by using this model [3] , and a model of a harmonic oscillator in an external 
gravitational field is considered and the developed formalism is applied to the study 
of thermal properties of noninteracting Bose and Fermi gases in harmonic traps [4-6]. 
It is necessary to note the generalization of this problem to the two-dimensional case, 
where an elastically bound particle in a space with a combined linear topological defect 
is studied in detail [7,8]. 

§ E-mail: ejafarov@physics.ab.az 



The Wigner d.f. of the relativistic finite- difference oscillator in an external field 2 



However, the solution to the problem under consideration in the configuration 
space is not enough and an important step here is to study the quantum problem with 
correspondence to its classical analogue. The phase-space approach allows us to answer 
this question and opens a lot of hidden features of the quantum system expressing it in 
the language of the classical approach. In order to see them, we just need to compute the 
quasi-probability function of the joint distribution of the momentum and position, where 
the Wigner distribution function is most known. Analysing the Wigner distribution 
function for the non-relativistic harmonic oscillator in an external field, we can observe 
that the applied linear field just shifts its stationary states to the negative position 
direction and does not have any influence to the values of the momentum. In this case, 
an important question about any role of the relativistic effect arises. The purpose of 
this paper is the phase-space study of the one-dimensional relativistic oscillator model 
in an external homogeneous field expressed through finite-difference Schrodinger-like 
equation [9-11]. 

Our paper has following structure. In section 2, we give brief information about 
the Wigner distribution function. The Section 3 is devoted to the non-relativistic linear 
oscillator in an external field, where we present its solutions both in configuration and 
phase spaces as well as the expression of the Wigner function of the thermodynamic 
equilibrium. Section 4 is devoted to the model of the finite-difference relativistic 
oscillator in an external field and in section 5, we present the explicit expressions of the 
Wigner distribution function of the stationary and thermodynamic equilibrium states 
for the relativistic linear oscillator in an external field. 



2. The Wigner quasi-probability distribution function 

The Wigner function W (p, x; t) [12] being as analogue of the classical distribution 

function in the phase space p (p, x) has a wide range of applications in the non-relativistic 

quantum mechanics [13-15]. The limit relation lim W(p, x; t) = p(p,x) between the 

h— »o 

Wigner function and the classical distribution function exists and it allows us to calculate 
quantum corrections to the known classical results employing an analytical form of the 
Wigner function. It is necessary to note that the first application of this function was 
the calculation of the quantum corrections to the classical distribution function of the 
equilibrium states of particle system in an arbitrary potential field. The Wigner d.f. 
is a function of momentum p and position x as well as a fucntion of the time t in a 
general case. One can obtain it from the wavefunctions of the quantum system under 
consideration both in position representation ip (x) and momentum representation (p) 
by using well-known transitions: 



W (p, x; t) 



2^/^^ + r' t ) e '"' V ( X "^'*) dl ' 



(2.1) 
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W (p, x; t) = i J 0* (p + ~p', tj e^fy (v - \p\ tj dp'. (2.2) 

-co 

The Wigner function (I2.1l) - (l2.2p satisfies the following equations: 



W(p,x;t)dp = \ip(x,t)\ 2 = W{x,t), (2.3) 
W(p,x;t)dx = \<j)(p,t)\ 2 = W(p,t). (2.4) 

Here, W(x,t) is the probability of the particle observation at point x at time t. 
Correspondingly, W(p, t) is the probability of the particle observation in the momentum 
space with the momentum value p at time t. 

Despite the fact that the Wigner function W (p, x; t) satisfies equations (12.31) and 
( 12. 4ft , one cannot consider it as the probability of the observation of the particle with 
the momentum value p at point x due to the fact that the function W (p, x; t) at some 
values of the p and x can become negative. 

With the help of the Wigner function one can find the average value of any physical 
parameter f(p,x) through the following formula: 

f — J f(p,x)W(p,x,t)dpdx. (2.5) 

An explicit expression of the Wigner function has been already obtained for a 
number of the non-relativistic quantum mechanics problems [13-17], whereas the phase 
space of the relativistic model of the linear oscillator described by the finite-difference 
equation [18] was considered in [19]. 

The purpose of this paper is to obtain the explicit expression of the Wigner function 
for the relativistic model of the linear oscillator in an external homogeneous field. 

3. The non-relativistic linear oscillator in an external field 

Let us first to consider in short the non-relativistic case. In the non-relativistic quantum 
mechanics, the Hamiltonian of the linear oscillator in a homogeneous external field 

TT0 h 2 d 2 mu 2 x 2 , n , 

has the following eigenfunctions [2]: 



^Nn^) = °Nn " H n ^(x + X ) 



a^L) .e-^aW (3 .2) 
h 



These wavefunctions correspond to the energy spectrum 



E 9 Nn = hnj(n+ ± -)-'^xl (3.3) 
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with Xq = g/muo 2 . Function (13. 2p satisfies the following orthonormalization relation: 

^Nn( x )^Nm( x ) dx = <W (3.4) 



From (13.41) it follows that 



W CN °-\^h) 



c Nn z , 

V2"n\ 

One can obtain the wavefunctions ip% n {x) (13.21) by the simple transition from the 
wave functions ip° Nn (x) of the non-relativistic oscillator at g — 0, i.e.: 

Vft„(*) = e-^r Nn ( x ) = ^ Nn ( x + x o)- 
In the momentum representation this transition is as follows: 

<t>Nn(P) = e ^<t>Nn(p)- 

After substitution of (13. 2p into (12.11) one can perform the integration and we will 
find the Wigner function of the non-relativistic linear oscillator stationary states in an 
external field (13. lft : 

W 9 Nn (p,x) = t}L e -(n^ f) Ln (2r? 2 + 2(£ + £o) 2 ) , (3.5) 



where ri = , p = and f = x x f^ are the dimensionless variables, £ Q = £n\/^r = 

1 vmuJh s V h ' ^ u u V h u V h 

is the dimensionless parameter and L n (x) is the Laguerre polynomial. 

By using equation (12.21) . one can obtain the following operator form of (13.51) : 

W Nn (P, x) = ^^H n (ri + l -d^j .H n (rj- e^^l . (3.6) 

It is necessary to note that the Wigner function (13.51) is normalized by the condition 

W 9 Nn (p,x)dpdx = 1, 

from where we obtain that 

^ r(fc + l/2)r(n-fc + l/2) 

h = ( } 

Another summation formula that we will obtain through the substitution of (13.21) 
to ([211D is: 

n 

2 k C k n ■(2k-l)\\- H 2n _ 2k (xV2) = 2 n H 2 n (x) , (3.8) 

k=0 
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Figure 1. The behaviour of the ground-state Wigner function of the non-relativistic 
linear oscillator in an external homogeneous field for values of the external field g = 
and 3 = 3 (m = u) = % = 1) 



where by definition 0!! = (—1)!! = 1 [20]. If we take into account the equality 
(2k — 1)!! = A^T (k + 1/2), then this formula can be written in an equivalent form 

as 

J2 2 2k C k n T (k + 1/2) H 2n _ 2k [xVfj = TRl (x) . (3.9) 

The Wigner function of the quantum system in the state of the thermodynamic 
equilibrium at temperature T is determined by the formula 

W 9 N ( P ,X) =Y, W NnW 9 Nn (p,x), (3.10) 



n=0 



where 



<n = ^7^ = 2sinh(^/2)e = w% w (3=l/kT. (3.11) 

It is possible to perform the summation f 1 3 . 1 j) for the case of the non-relativistic 
linear oscillator in an external field, where we will obtain an explicit expression of the 
equilibrium Wigner function as follows: 

W£(P, x) = tanh ( ^ /2) exp [- (r? 2 + (£ + ^ ) 2 ) tanh(^/2)] . (3.12) 

The behavour of the ground state Wigner function Wjy (p, x) determined by 
equation (13. 5p is presented in Figure [H On can see that the applied external field 
in the non-relativistic case does not influence the joint distribution form just shifting it 
along the negative values of the position. 
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4. The relativistic linear oscillator in an external field 

The relativistic model of the linear oscillator in an external homogeneous field is 
described by the finite-difference Hamiltonian [9-11] 

2 

H 9 (x) = mc 2 cosh iXd x + -^-x(x + iX)e iXdx + gx, (4.1) 

where A = h/mc is the Compton wavelength of the particle with mass m and parameter 
c is the speed of light. 

Hence, one needs to mention that the dynamical symmetry group and Barut- 
Girardello coherent states were constructed for model (14. ip and the bilinear generating 
function for the Meixner-Pollaczek polynomials is obtained in [9] and the generalized 
coherent states for the model of the relativistic linear oscillator were constructed in [11]. 

The eigenf unctions corresponding to the Hamiltonian H 9 (x) ( 14. ip have the following 
form: 

V*(*) = 4 (^) " * T(u + ix /X) P v n (|; y) e (*-*/*W* (4.2) 
where the normalization constant c 9 




and P^{x; ip) is the Meixner-Pollaczek polynomial. 

The wavefunctions (14. 2 j) are orthonormalized by the following condition: 

i^'(x)'^~l n (x)dx = 5 nm . 

The energy spectrum of the relativistic oscillator in an external field corresponding 
to the Hamiltonian H 9 (x) (14. ip and the wavefunctions (14.21) has the following form: 

E 9 n = hw5 (n + u) , 71 = 0,1,2,3,..., (4.4) 

where the spectrum (I4.4p becomes discrete under the condition \g\ < mcuo and the angle 
< (p < tc is determined through the relations cos<£> = g/mcu and 5 = sin (p. 

One can obtain the wavefunctions <fi 9 n (p) of the quantum system in the momentum 
representation from the wavefunctions ip^x) (14.21) through the relativistic Fourier 
transform: 



{p) = -j= j C(p,^ 9 n (x)dx, (4.5) 



where the function 



\ ix/X 

e(p,x)=(^±-) =e^ (4.6) 



mc 
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is the relativistic plane wave [18,21] and \ = ln((pQ+p)/mc) is the rapidity. As a result, 
we find that the wavefunctions of the momentum representation are expressed through 
the generalized Laguerre polynomial: 



(p)=c> n Ce^Ll»- 1 (26C), 



(4.7) 



i n (26 f 



nl 



mcT(n + 2v) 



c(po + p) 



mc 
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5. Wigner representation of the relativistic oscillator in an external filed 

Let us now construct the Wigner representation of the relativistic oscillator in an 
external homogeneouf field. First, we will consider Wigner function of the oscillator 
stationary states and we will use for our calculations the following definition of the 
Wigner function for the relativistic case [19,22]: 



W»(p, x) = JL J 0f (x + ~X') <t>l (x - \x^j e-^ 7 W, (5.1) 

— oo 

where <p 9 n (p) = 

Substituting (14.71) into (15.1 ft and taking into account the relation 

±Sa 2ix\' I-,/ 2ixx' 

e ± ~ ax e = e x e 

we will obtain the following finite-difference expression for the Wigner function of the 
stationary oscillator states in an external field: 

W'fax) = ^/n^me^L 2 :- 1 (25Ce^) W° (p,x) . (5.2) 

Now, we can obtain the explicit expression of the Wigner function of the relativistic 
oscillator ground state in an external field as follows: 



oo 

\2v 



With the help of the following integral relation [20] 

a/2 



x a - l e~ px - qlx dx = 2 m K a {2y/pq), Rep > 0, Reg > 0, 



(5.4) 



one can express ( 15.31) through the Macdonald functions: 



w ° (P) x) = lww/ 2 ^ )x/XK2lx/x (2 ° ' (5 ' 5) 
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In a similar way, we obtain the following expression for the Wigner function of the 
excited states: 

W 9( v x ) - 2(2<5Q 2 ^! (2y _, w> 

* £ ( r r 1 ) ( 1 ) ^^^). 

where (a) = a (a + 1) • • • (a + n — 1) = r (a + n) /r (a) is the Pochhammer symbol, 
r(a+i) 



a 
n 

of the generalized Laguerre polynomials: 



r( a +i-n) are ^ ne binomial coefficients and here we also used the explicit form 



fc=0 ' \ 



n + a ^ k 
n — k 



The limit relation, reducing the (15.61) to the (13.51) is presented in detail in the 
appendix. 

As in the case of the relativistic oscillator in absence of the external field [19], we can 
determine the equilibrium Wigner function for the relativistic oscillator in an external 
field by formula (13.101) with the coefficients 

w g n = e -PEi/ Z 9{(3), (5.7) 

where E% is the energy spectrum of the relativistic oscillator (I4.4p and the partition 
function has the following form: 

Z»(/3)=;>-« ^ (5.8) 

n=0 

Let us express the equilibrium Wigner function by the following series of the 
integrals: 



1 oo 

W " ( P > X ) = ~7q(~R\ E e ^ KW n(P, *) (5-9) 

oo 

1 _ p -PhuS °° \ J3hjj&n [■ 

■ O 2 *) E r(n + 2v) J e a ^ y+a ^- y Ll^ (25^) (2^) e**"*dy. 



nh 

n=0 



-oo 



One can now use the following bilinear generating function for the generalized 
Laguerre polynomials [23] 

E r( n + a + i) L - (a;)L " (z/) = "T^- 7 " yY^~r ) exp ' N < 1(5 ' 10) 
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Figure 2. The behaviour of the ground-state Wigner function of the relativistic linear 
oscillator in an external homogeneous field for of speed of light c = 4 and values of the 
external field g = 0, 1, 2 and 3 (m = u> = H = 1) 



and changing the integration and summation for the function W 9 (p, x) (15.91) we obtain 
the following expression for the equilibrium Wigner function: 

W(p,x) = ^e-^/V-^/^ (~r^) ^C^ 2 + 5 2 coth 2 ^,(5 

where 4> = arg (ip + 5coth ^y^), p = cos<p as well as I a (x) is the modified Bessel 
function of the first kind. To derive (15. Ill) , we have also used the formula (13. 6ft . 

It is obvious that the obtained expressions of the Wigner function for the relativistic 
oscillator in an external field in case of the g = (i.e., 5=1 and p = 0) go to the 
corresponding expressions of the relativistic oscillator without external field [19]. 

The behaviour of the ground-state Wigner function of the relativistic linear 
oscillator in an external homogeneous field Wq(p,x) expressed through equation (15.51) 
is presented in Figure [2j One can mention here that the limit c — > oo recovers the non- 
relativistic case and behaviour of ( 15. 5ft will be identical to the non-relativistic ground 
state Wigner function Wjy (p, x) determined by equation (13. 5p . which is presented in 
Figure [U Then, any finite value of speed of light c should show for us the contribution 
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of the relativistic effects. Therefore, we depicted the joint distribution for the value of 
speed of light c = 4 (m = u = h = 1) and see that unlike the non-relativistic case, more 
strong field to be compared with the speed of light cardinally changes its Gaussian-like 
form as well as shifts this distribution along the negative values of the position. One can 
see that the distribution in the absence of the field is determined for both positive and 
negative values of the momentum whereas a more strong field shifts it to the positive 
momentum values. 

We do not present here graphical distribution of the Wigner function for the 
thermodynamic equilibrium state and the excited states of the relativistic oscillator 
under influence of the external field. Only we need to note here that the behaviour of 
the relativistic oscillator excited states in the phase space is similar to the ground-state 
Wigner function behaviour. Same picture can be observed for the Wigner function for 
the thermodynamic equilibrium states, which is a good approximation of the ground 
state joint distribution for the room temperature. 

Appendix A. 



In this appendix, we compute the non-relativistic limit of the Wigner functions of the 
relativistic linear oscillator in an external field. First, we have to take into account that 



p = COS If = 

v 7 / 1 



, 5 = siny? = y/l - ^, x/X = ( = yf\L (r\ + vV + ?? 2 ) , 



mc 2 

Then, at fi — > oo we will have the following limit relations for the elementary 
functions: 

lim 2(f)x/X = —7^, 

IX— »oo 

lim (2f - tt) x/X = -2ff , (A.l) 

fl— >oo 

as well as for the special functions: 

h,-! (2ioC) = _i_e 2 ""-«+ 2 ^"+(»'^)" 2 -'«8 j (A.2) 

V47T/iO 

2 

7T ^9,,h—0 mh^—hr^ — L. 



K 2ix/X (2k) = yj e -^VPM* 
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where: 

To obtain the formulae ( 1A.2I) we used the following representations for the modified 
Bessel, the Macdonald and the gamma functions [24]: 



e 



^Jp 2 +x 2 —p sinh 1 £ 



2nl p (x)^ , x^p^oo 



x 2 



K ip ^) = yj\- • e-V^-P^ 1 S, x > p > 0, p _> oo, 

as well as one needed to take into account the following limit relation between the 
Laguerre and the Hermite polynomials [25]: 

(-l) n 

lim zz-tL^" 1 (2u + 2y/ux) = ^-/-H n (x) . 

By using the equations ( 1A.1I) and ( 1A.2I) . one can show that at c — > oo, expressions 
(15. 2p and (15.111) coincide with the expressions (13.61) and (I3.12p . respectively. 
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